apparatus and method for detecting flinch 
Response to Nociceptive Agents 



BACKGROUND OF THE INVENTION 

1. Field of the Invention 

^ The present invention concerns detection of response to nociceptive pain. 

This invention particularly concerns an automated model for screening anti- 
nociceptive agents by detecting and measuring a flinch response in an animal 
whose extremity has been subjected to an irritant such as formalin. 

2. Description of the Prior Art 

1 0 The escape response or agitation evoked by a transient, strong stimulus 

attests to there being a close relationship between stimulus intensity, peripheral 
afferent discharge and magnitude of the pain state as defined by response latency 
and magnitude. There are situations, however, in which the magnitude of the 
response to pain may exceed what would normally be anticipated given the 

1 5 magnitude of the physical stimulus and the afferent traffic generated by that 

stimulus (31, 45, 47). These situations are loosely considered as reflecting a state of 
hyperalgesia, possibly arising from sensitization of the peripheral terminal aid/or a 
central facilitation. 

Several preclinical models have been developed that may reflect the 
20 significance played by such facilitation on behavior. The common characteristic 
found in these models is the injury that is induced and its causing of the sensory 
axon to produce a persistent discharge. A frequently used method of producing 
injury in the rat is the subcutaneous injection of a small volume of irritant such as 
formalin into its hind paw. Typically, after the formalin injection, the rat displays a 
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biphasic (Phase I and Phase II) incidence of flinching (rapid paw shaking) and 
licking of the injected paw (18, 42, 43). The behavioral syndrome produced by the 
injection of formalin into the paw has been widely used to define the pharmacology 
of systems that regulate facilitated processing. The "formalin test" has evolved into 
5 a widely used tool in the screening of analgesic and anti-hyperalgesic drugs (45). 

An important limitation of this behavioral model is its labor-intensive nature 
regarding data collection and the time required to train observers in its reliable 
implementation. Several automated systems have been proposed to facilitate data 
collection. One approach has been to employ strain gauges to measure mass 
1 0 movements of a rat in a confined cylindrical cage (2 1 ). A second model involves a 
video system that employs a pattern recognition algorithm (22). Though of merit, 
these approaches are limited in that each only indirectly measure movement of the 
injected paw. Therefore, there is a need for an approach that directly measures 
movement. 

1 5 This invention addresses this limitation by presenting a system that assesses 

only the movement of the injected paw. As described below, this approach involves 
placing a metal band on the injected paw and detecting the movement of that band 
with a localized low strength sinusoidal electromagnetic field. 

SUMMARY OF THE INVENTION 

20 

One object of the present invention contemplates the capability of 
electronically detecting a response to chemically induced nociceptive pain directly. 

Another object of the present invention contemplates a method for rapid 
screenning for nociceptive agents, using an automated model. 

To achieve these objectives, the most preferred embodiment of this invention 
25 is an automated flinch-detection apparatus for measuring spatial displacement of an 
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animal's paw injected with irritant. This apparatus comprises an electromagnetic 
detecting assembly having (1) a transmitting occilator for generating electrical 
current (2) an electromagnetic transmitter coil coupled to the occilator for 
generating an electromagnetic field (3) an electromagnetic receiving coil placed in 
5 axial plane directly below the transmitter coil (4) a first, receiving amplifier 

connected to the receiving coil (5) an amplitude detector connected to the receiving 
amplifier (6) a second amplifier connected to the amplitude detector (7) a metal 
object attached to the animal's paw (8) and a cylindrical observation chamber of a 
diameter not greater than the diameter of the generated magnetic field. 

1 0 The chamber is placed directly over the receiving/transmitting coil assembly, 

wherein the current generated by the transmitting occilator circulates in the 
transmitter coil, creating an electromagnetic field that penetrates the metal object, 
creating eddy currents perturbing the electromagnetic field. The fluctuating 
perturbations are picked up by the receiving coil, amplified by the receiving 

1 5 amplifier, detected by the amplitude detector and further amplified, filtered and 
digitized. 

The current passing through the transmitter coil typically generates an 
electromagnetic field in the 6 to 8 kilohertz range with a signal strength on the 
order of 5 to 8 milliwatts. The metal object attached to the animal's paw is a small 

20 metal annular collar that fits around the paw. The object can also be a small metal 
"C" collar in incomplete annular form. The metal object can be made of, or contain 
a ferrous metal. The metal object can also be non-ferrous. The observation 
chamber is a transparent cylindrical container, insuring that the animal will remain 
inside the boundaries of the electromagnetic field generated by the coil assembly. 

25 The observation chamber may also individual compartments, permitting testing of a 
plurality of animals The detection assembly below the observation chamber, in this 
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case, would have multiple independent detection units. The observation chamber 
typically is constructed of any rigid transparent plastic such as, but not limited to, 
Plexiglas. 

Another most preferred embodiment contemplated is a method for ^^^u™g 
5 a flinch response by an animal whose paw has been subjected to an irritant. This 
embodiment comprises attaching a metal object to the animal's paw and placing the 
animal in an observation chamber, as described above, situated directly over a 
detection assembly having a transmitting occilator for generating electrical current 
and an electromagnetic transmitter coil coupled to the occilator for generating an 
10 electromagnetic field. An electrical current is thus generated and received by an 
electromagnetic receiving coil, which is connected to a receiving amplifier that 
amplifies the received generated electrical current. An amplitude detector receives 
the amplified electrical current, which, in turn, is amplified by a second amplifier. 

Thus, the current generated by the transmitting occilator circulates in the 
1 5 transmitter coil, creating an electromagnetic field that penetrates the metal object 
attached to the animal's paw, creating fluctuating eddy currents perturbing the 
electromagnetic field. The fluctuating perturbations are picked up by the receiving 
coil, amplified by the receiving amplifier, and detected by the amplitude detector. 
The pertubations are further amplified, filtered and digitized to produce a measured 
20 response to the applied irritant. 

These objects and other aspects and attributes of the present invention will 
become increasingly clear upon reference to the following drawings and 
accompanying specification. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1. Top: Block diagram of the detection unit and its functional 
elements. Bottom: Raw signal tracing over 2 minutes of an animal displaying flinch 
behavior with arrows added where the processing algorithm detected flinches. 

5 Figure 2. Photograph showing the four detection units in which rats are 

placed for flinch test and the analog input/signal processing module. The computer 
used in flinch detection, display and analysis is located below the testing table. 

Figure 3. T^: Displays the modeled transmitter (inner) and receiver 
(outer) coil configuration and modeled position of the rat paw collar. 

1 0 Bottom . The graph displays the variation of the magnitude of the magnetic 

field when in phase with the transmitter current across the bore of the coils using a 
fixed paw collar configuration. The Y axis presents the dimensionless ratio of "B 
(magnetic flux density at current location along the radius from coil center)" 
divided by "B center (magnetic flux density generated at the coil center)" The X 

1 5 axis shows the ratio of flux measured at the specified location along the radius 
divided by the radius of the receiver (outer) coil. Where the X- axis ratio is 1 .0 
the measurement is over the receiver coil. In the graph displayed, the transmitter 
coil is configured to be half the radius of the receiver coil. 

Figure 4. Representation of current flow occurring with open and closed 
20 collar configurations. 

Figure 5. Photograph displays the screen present during concurrent data 
collection in 4 rats. The view, for each animal, from left to right is: 1) 
Study/Animal. 2) The small green virtual "On/Off' switch that initiates and ceases 
data collection. 3) Smoothed range waveform signal of previous 2 seconds with 
25 indicators for flinch activity and flinch count for that interval. 4) Minute by minute 
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flinch count sums over completed portion of the study for each animaJ. 

Figure 6. Drawing indicates the paw band and the typical placement on the 
left hind paw. 

Figure 7. Toe: Average (+SEM) time effect curve for flinching behavior 
5 expressed as mean flinches for 100 rats. Middle ; Frequency of flinching counts for 
Phases I and II displayed as a standardized distribution (Z-Scores). Bottom : The 
mean cumulative flinches (+SEM) observed during Phase I (0-9 minutes). Phase II 
(10-40 minutes). Phase IIA (10-40 minutes) and Phase EDB (41-60 minutes). 

Figure 8. Control Charts. The Xbar Line Chart provides information about 
1 0 variation between the subgroup means over time by displaying individual group 
mean flinch counts plotted along its time (Date) axis. The S Line Chart plots 
subgroup standard deviation along the time axis providing information about 
variation within the subgroups. Calculated upper and lower control limit (U/LCL) 
and average (Center) lines have been added. 

1 5 Figure 9. Time effect plot comparing flinches per minute as determined by a 

trained human observer and by computer (Left) and a scattergram (Right) showing 
linear regression between the concurrent human and computer flinch detection of 
six rats. 

Figure 10. Time effect curves for vehicle, intrathecal (top left) and 
20 intraperitoneal (top right) morphine on formalin induced flinching in the rat. 
Cumulative flinches, by dose, are presented by phase for intrathecal morphine 
(bottom left) and for intraperitoneal morphine (bottom right). 

* p<0.05, 1 Way ANOVA across dose. 

Figure IL Time effect curves for vehicle, intrathecal (top left) and 
25 intraperitoneal (top right) MK801 on formalin induced flinching in the rat. 



Cumulative flinches, by dose, are presented by phase for intrathecal MK801 
(bottom left) and for intraperitoneal MK801 (bottom right). 

* p<0.05, 1 Way ANOVA across dose. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

5 Description of the System. 

The detection device consists of a pair of electromagnetic coils as depicted 
in Figure I, one serving as a transmitter 2 and the other as receiver 4. When 
current passes through the transmitter coil 2, an electromagnetic field in the 6 to 8 
kilohertz range is generated with a signal strength on the order of 5 to 8 milliwatts. 

10 Eddie currents are set up in ferrous or non-ferrous metals within this field and it is 
movement of the metal in the field that is detected. Location, size, electrical 
conductivity and magnetic permeability are factors that determine how well 
movement of the metal is detected. The receiver coil output is amplified 10-14, 
filtered and digitized for analysis. The spatial displacement of the injected paw 

15 within the electromagnetic field is detected using a small metal collar placed on the 
paw, as seen in Figure 6.^ During testing, the animal is placed into a cylindrical 
Plexiglas container 6, Figure 1, (15 cm diameter x 30.5 cm high) mounted above 
the transmitter/receiver coils assembly which is contained within a plastic 
enclosure. The cylinder 6 insures that the animal (rat) will remain inside the 

20 electromagnetic field generated by the coil 2 without the added stress of being in 
restraint. As shown in Figure 2, the system includes independent detection units 
and permits concurrent testing of 4 rats in this case. 

Modeling of physical system parameters 

Several parameters of the electromagnetic field and paw band interaction 
25 were initially modeled to understand operating characteristics of the flinch 



detection system with the aim of defining the contributions of 1) collar 
permeability; 2) a complete or "0"-shaped collar configiiration versus a partial or 
"C"-shaped collar configuration; and 3) relationship of collar mass and sinusoidal 
fi*equency. 

5 Modeling assumptions . 

The following modeling assumptions were made in order to define the 
effects of the several variables on system performance. 

1) While the collar may be variably located and oriented during normal 
study conditions, for the puipose of this evaluation, the collar was assumed to be 

1 0 centered with respect to the transmitter and receiver coils, thus simplifying the 
model to an axisymmetric representation of the collar and coils. The coil motion 
was limited to moving along the coil axis in a vertical direction (see Figure 3), 

2) The collar was considered to be permeable and the representative values 
were assumed to range JBrom 100 to 1000 (low-grade ferrous) or 100,000 

1 5 (permaloy). It was assumed that the system contained only materials that could be 
represented by a single permeability value. 

3) Given the dimensions of the coils, the magnitude and the fi-equency of 
the voltage (e.g., wavelength relative to coil diameter), the capacitive effect was 
considered to make a negligible contribution. 

20 4) Motion of the collar was assumed to be limited to height range of 

between 0.35 and 0.45 inches above and perpendicular to the coils; and periodic, 
with a period of motion (between 1 and 20 Hz., see below). It is important to note 
that the field sinusoidal fi-equency was chosen to be large compared to the target 
motion (6-9 kilohertz). The extreme positions of 0.35 and 0.45 inches were 

25 selected for this modeling to correspond to the elevation extremes of the raf s paw 



above the coils. Even though the collar is changing its position, this assumption 
permits the motion of the collar to be treated as a steady state condition. While the 
direction of motion and coil dimensions will affect the receiver output, the effects 
due to the collar material, collar thickness and the open vs. intact configuration of 
5 the collar can be identified by using a simplified model. 

5) For modeling purposes, antennae coils were considered to be of a fixed 
loop design and constructed of small strands to permit minimal eddy current. 

6) The properties of the paw collar (Figure 4) were varied in terms of its 
configuration ("C" or "O"), thickness and material permeability. Collar thickness 

10 was of interest as it is related to collar weight and the ability of the animal to move 
freely. Magnetic permeability also was considered to be a significant parameter as 
it governed whether or not collars were to be made of specialized materials. These 
factors affect the distribution of the current generated in the collar and its magnetic 
flux. Collar configuration was important as it related to construction and mounting 

15 and distribution of the currents 40 generated within the collar by the transmitter 
coil. A continuous circle 42 tends to collect maximum current flow along its outer 
diameter while an incomplete circle 44 (the arc) current distribution loops around 
the entire surface perimeter as shown in Figure 4. 

6) The electromotive force (EMF), in volts, appearing across the receiver 
20 coil is the time derivative of the magnetic flux (<|)), in Webers, being linked with 
the receiver coil of N turns. Since the supply signal is sinusoidal at a frequency (f), 
in Hertz, the EMF across the receiver coil was computed by the formula 2%(N^, 
where f = 8KHz. Using the assumption that the paw motion can be treated as a 
steady state problem, the maximum and minimum collar elevations were examined. 
25 Electromagnetic analysis, employing the sinusoidal assumption, was performed for 
each position of the collar. Two solutions were produced: 1) an imaginary solution 
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which represented the electromagnetic field 90° out of phase with the input signal; 
and 2) the real solution where the electromagnetic field was in phase with the input 
signal. The peak EMF was computed by combining the EMF for each solution 
using the square root of the sum of the squares. The effect of collar motion is 
5 represented by calculating the differences of EMF at the two extremes of collar 
position. Calculations were carried out using the commercially available software 
ANSYS, employing the finite element method to represent the electromagnetic 
field. Guidelines for analysis of sinusoidal varying fields are provided with the 
software. 

10 Modeling results . 

The modeling reflects the analysis of four combinations from two collar 
thicknesses, 0.05 inch and 0.005 inch, and two collar configurations, open and 
closed circumference. These four cases were evaluated using a variety of magnetic 
permeability's of the collar (100-100,000) and are summarized in Table 1. The 
i 5 results are normalized by a single case. The values shown represent the relative 
behavior of the system with respect to alterations in collar design. It is understood 
that EMF can be also altered by the number of turns in the receiver coil or the 
transmitter coil and larger paw motions would result in larger values, but the trend, 
shown by these results, would not be altered. 

20 As noted above, the data presented in Table 1 reflects the analysis relevant 

to the paw collar being in the center of the antenna coil. The relative contribution 
of the several variables would, however, be the same if assessed in any part of the 
field. To determine the EMF generated across the field, the analysis was carried 
out for a standard coil across it's radius where the transmitter coil is constructed to 

^5 be half the radius of the receiver coil. A typical variation of the normalized field 
across the bore of the transmitter and antenna is shown in Figure 3. Values larger 
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than unity correspond to the magnetic flux density (B) being larger than the 
magnetic flux density at the centerline of the coils. This variation is normally 
observed for such coil configurations. Larger flux densities lead to larger values of 
EMF since the magnetic flux {^) is the area integral of the flux density over the 
bore of the antenna. Collar motion nearer the transmitter radius thus produces 
larger EMF. Outside the receiver coil, the field is observed to decay rapidly to a 
level at which a response would not be recorded regardless of the collar design. 

Table 1. Summary of collar modeling results for voltage signals generated by 
movement of paw collar in electromagnetic field generated by a surrogate flinch 
detection device. All EMF receiver values are normalized to Case 2, Pemieability of 
1,000. 



Permeability 


Case 1 
(.05", O-shape) 


Case 2 
(.05", C-shape) 


Case 3 
(.005", O-shape) 


Case 4 
(.005", C-shape) 


100 


0.07 


0.09 


0.06 


0.06 


150 


0.11 


0.13 


0.09 


0.09 


175 


0.79 


0.90 


0.64 


0.64 


250 


0.86 


0.95 


0.64 


0.64 


500 


0.88 


0.98 


0.67 


0.71 


1,000 


0.95 


1.00 


0.69 


0.76 


100,000 


1.05 


1.05 


0.76 


0.76 



Modeling Conclusions . 

This analysis provides several conclusions that were of concern during 
1 5 system optimization. 

L The most significant factor was metal permeability. Typical permeability for 
standard iron is around 1,000. It should also be noted that, if collar cross sectional area is 
reduced, permeability would be further reduced. The significant drop in the values at a 
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permeability level of 150 and corresponds to the physical conditions of a collar where 
current cannot be concentrated at its surface. Additional sensitivity could be derived from 
the use of materials with high permeability and/or with increased thickness. While Cases 
3 and 4 show a 20% to 30% decrease in the relative performance, it must be remembered 
5 that the mass of the thinner collar is only 10% of the thicker collar. 

2. The analysis confirms that use of an open collar does not result in major 
degradation-of system performance. 

3. The use of higher permeability materials (>1000) reduces concern with 
collar geometry. 

^ ^ 4. The relationship between frequency of sinusoidal current flow and generated 

EMF is linear, so increasing energy will increase signal and system sensitivity. 

5. Modeling the cross sectional field profile makes clear that the field 
intensity peaks over the antennae coil. The area to which liie animal is restrained 
should thus be centered over this radius. Importantly, profiling emphasizes the 

1 5 rapid tailing of the flux density as one extends beyond the radius of the receiver 
coil. This serves to limit the external influences of the field. This, plus the low 
frequency /energy signal produced by this system reduces the likelihood of EMF 
interference with other electronic devices. 

6. For any given transmitter frequency and energy level, the flux density will be 
20 increased by reducing the relative diameter of the transmitter and receiver loops. Thus, 

the magnitude of the EMF generated for a given paw movement and collar parameter can 
be increased if enhanced sensitivity is required (e.g., as with a smaller collar or smaller 
paw displacement). 
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Signal analysis 

Signal conditioning . 

The analogue signal obtained (Figure 1) from the receiving coil 2 is filtered at 3 
hertz and amplified 10-14 before being digitized at a sampling rate oflOOO Hertz using 12 
bit resolution. Beside the flinch behavior that occurs during the testing process, any 
general displacement of the limb, such as during ambulating and/or grooming activity, 
also generates signals that need to be addressed To minimize the interference caused by 
non-flinch movement, the digitized signal is subjected to real time analysis using a 
software (LabView) algorithm to pick out a flinch from other paw movement. 

Signal detection . 

Two approaches were considered in selection of triggering algorithm. The first 
employed a "power/frequency" analysis, in which spectral analysis of the limb's flinching 
movement revealed a characteristic frequency at about 8 Hertz. The peak in the power 
spectrum was found, however, to be very broad, with significant components observed as 
low asl Hertz and as high as 20 Hertz. The breadth of the flinch bandwidth caused the 
algorithm to be judged not sufficiently discriminatory at differentiating flinching behavior 
from regular paw movements. 

The second approach was a "zero crossing interval/peak heighf analysis based on 
estimating the amplitude of the signal over a sliding time window. The range of the 
voltages (maximum-minimum differences) over a moving 128 millisecond interval was 
calculated to produce a continuous output waveform. This secondary or range waveform 
contained jagged peaks that correlated well with the flinch transients found in the 
acquired waveform. It was then smoothed, using a linear convolution filter (a 128 
millisecond non-weighted moving average filter), and the smoothed range waveform was 
examined in real time by a peak detection algorithm set to pick out spikes of -300 
milliseconds duration and amplitudes of > 0.5 volts. Each peak thus detected was counted 
as a paw flinch. This algorithm was observed to produce flinch reports that correlated well 



13 



with data obtained from animals in which concurrent scoring by trained observers were 
obtained (see below). 

Data collectioii 

Signal events that meet the criteria of a flinch, as explained above, are 
5 captured and summed by time interval (normally set at one-minute increments) over 
the course of the study. An example of a typical "real time" signal showing flinch 
behavior is displayed in Figure 1 (bottom). Signals are collected and stored in 
master raw data spreadsheets by animal according to its Study/Animal identifier. 
The software (Labview and VisualBasic) presents the data for 4 animals 

10 simultaneously on the screen (Figure 5) and includes for each animal: a) the 
Study/Animal identifiers including animal numbers, treatment codes, dates and 
other information relevant to the study; b) a window displaying the previous 2 
seconds of digitized signal with markers indicating any flinch detection activity and 
a count of the flinches within that window; and, c) a line graph of each animal's 

1 5 flinch count by sampling the interval (1 minute) from initiation of testing. A typical 
view of the computer screen during testing of 4 animals is shown in Figure 5. 

Preparation of the animal 

A soft metal band (10 millimeters wide by 27 millimeters long, shaped into a 
"C", and weighing approximately 0.5 gram) is placed on the hind paw of the animal 

20 being tested. The open part of the "C" is positioned at the top of the paw with the 
arms of the "C" gently compressed to form a bracelet around the paw (see Figure 
6). Retention of the band is enhanced by applying a small amount of adhesive 
(Cyanoacrylate, Elmers, Columbus, OH). The size and weight of the band is 
sufficiently small so as not to hinder the animaFs normal movement. Animals are 

25 allowed to accommodate in individual Plexiglas chambers for 1 hour before being 
moved to a test chamber. Just prior to the animal being placed into the test 
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chamber, it is briefly restrained in a cloth towel, and an irritant (typically 5% 
formalin, in volumes of 50 is injected into the dorsum of the banded paw. Data 
collection is initiated after the animal is placed inside the test chamber. 

Drug delivery 

5 To examine the effects of a drug on flinching behavior, rats received 

intraperitoneal or intrathecal injections of the drug. Intraperitoneal injections were 
delivered in volumes of 0.5 ml/kg. Intrathecal injections were done in rats that had 
been previously implanted with chronic intrathecal catheters (see below) and using 
drug volumes of 10 |il followed by a 10 nl flush using the vehicle. 

10 Drugs/Chemicals 

The drugs delivered intrathecally and systemically were morphine sulfate 
(Malinkrodt) and MK801 (RBI) dissolved in physiological saline (0.9% NaCl w/v). 
The formalin solutions were prepared by diluting formaldehyde (Fisher Scientific 
Fair Lawn, NJ, Formalde-Fresh 20%) with physiological saline (0.9% NaCl w/v). 
1 5 Solutions were prepared fresh daily. 

Data analysis 

Primary data management. 

Each animal's flinch count value over an interval of time (usually 1 minute), 
for the duration of the study (usually 60 minutes), forms the data set used in all 
20 subsequent analyses. These data are averaged to characterize similar groups of 
animals by mean and standard deviation. For graphic display of data, flinches per 
minute are normally presented; but the mean of the flinches observed over 5 minute 
time periods (12 per hour) have been used to clarify the graph when a number of 
groups are being displayed for comparison of time effect. For statistical 



comparison, the total number of flinches observed during any selected time period 
(Phase I, n, DA, IIB) is calculated by accumulating each individual animal's flinches 
over that time period and averaging the group. Spreadsheet software (EXCEL), 
used in conjunction v^ith the data collection software, is constructed so as to allow 
5 the user to select up to 4 phases within the study's time frame for making phase- 
related analysis, including phase-cumulative flinch averages, standard deviations, 
and standard error values, and calculates Per Cent Maximum Possible Effect 
(%MPE), group means, standard deviations and standard error values when a group 
of control animals is available. Dose-response curves are based on the calculated 

10 %MPE. The flinch time course has frequently been divided into two principal 

phases: Phase I (0-9 min) and Phase n (10-60 min). In this facility it been observed 
that different drug effects may arise between early and late Phase 11, so, additional 
analysis occur at Phase HA (10-40 minutes) and Phase IIB (41-60 minutes) (23). 
For statistical comparisons, ANOVA is normally carried out between treatment 

1 5 groups for each of the 4 phases, with statistically significance-triggering post hoc 
test comparisons. Control groups normally consisted of a minimum of 5 rats run in 
temporal proximity to the drug groups. 

Dose response analysis . 

Dose response slopes and ED50 (effective dose in 50% of animals) values 
20 with 95% confidence intervals are calculated using a least squares linear regression 
analysis carried out on per cent inhibition of the cumulative flinching number for 
each phase (37). The per cent inhibition is computed by dividing the response of 
each animal in a given group by the respective vehicle control group and 
multiplying by 100. 

25 Statistical power . 

The minimum difference between the flinch count number measured in a test 



group and that measured in a control group, and determined to be statistically 
different, was calculated. The calculations were based on data obtained using 
independent control groups, while assuming different power estimates in the 0.7 to 
0.9 range, fixed group sizes (N = 8), and where power is defined as 1 - the 
5 probability of a Type 11 error being committed. Data were calculated as described 
elsewhere (52). 

System measurement stability . 

A program was constructed for estimating measurement expectations from the 
paw flinch system under control conditions, and for documenting long term process 

0 reliability. Data was collected at regular intervals over a period of 7 weeks for the 
purpose of preparing control charts to answer the foUov^ng quality control issues: a) 
what is the statistical distribution upon which we can define the process as being "in 
control," in other words stable and distributed normally; and b) does the process show 
itself to be "capable", does it meet assigned specification limits. Control and 

5 specification parameters (Upper and Lower Control Limits) were defined at study outset 
to be 3 standard deviations of the generated flinch count distribution for Phases I and n. 
Analysis was done using the Quality Control Tools found in Statview v.5.0 (SAS Institute 
Inc. SAS Campus Drive, Cary, NC 27513). 

The animals were acclimated and then treated only with formalin injection just 

0 before the start of the 60 minute data collection period. For each phase, the following 

values were calculated: a) the Xbar chart Center Line Value, an estimate of the process 

average and computed by averaging subgroup averages; b) the Xbar chart Control Limits, 

which identify the distribution range within which the process can be considered in 

control and computed as 3 times the estimate of the process standard deviation divided 

5 by the square root of the subgroup size (n = 4); c) the S chart Center Line Value, an 

estimate of the process variation computed by multiplying an unbiasing constant based 

upon subgroup size (n = 4) by the process standard deviation estimate; d) the S chart 

Control Limits, which identify the distribution range within which process variation can 
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be considered in control and computed as 3 times the estimate of the process standard 
deviation times an unbiasing constant based upon the subgroup size (n = 4); e) the % 
Upper and Lower Specification Limits that indicate the number of times the measured 
values exceeded either the upper or lower control limits and displayed as a % of the total 
number of measured values; and f) the Capability Index (Cp), an indication of whether 
the process is in control and can stay within the limits specified as necessary for test 
relevance, and calculated by dividing the Upper and Lower Specification difference by 6 
times the process standard deviation estimate. Upper and Lower Specification limits were 
calculated by phase as the mean of all tested animals ± 3 times the standard deviation of 
the tested population, where, if the Lower Specification limit was less than 0 that limit 
became 0. A Cp value of greater than L33 means that a very small number (6 out of 
100,000) of the formalin tests do not fall within the limits considered relevant to the 
normal testing process (e.g., Cp<LO, process not capable; Cp=LO, process marginally 
capable; Cp>LO process is capable). 

Animals 

Rats were male or female Sprague Dawley, Holtzman (Indianapohs, IN). Unless 
otherwise stated, the typical weight and age of these animals was 275-300 grams and 100- 
120 days. All animals were given a minimum acclimation period of 3 days on site prior to 
being entered into a study. For those studies requiring drug to be delivered spinally each 
rat was implanted with a chronic intrathecal catheter. The catheter, constructed of 
polyethylene tubing (PEIO), was implanted under halothane anesthesia by inserting it 
intrathecally though the cisternal membrane and passing it 8.5 cm. caudally to the rostral 
edge of the lumbar enlargement. It was then externalized percutaneously at the top of the 
head for access during drug delivery (49). Animals receiving intrathecal catheter implants 
were allowed 4 to 5 days of recovery before testing. No animal showing neural deficit or 
behavior abnormality was used in the studies. All animals were euthanized immediately 
after completion of testing. 
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study objectives, 

1 . Testing for "time of day'' and test chamber effects on flinch behavior. 
Adult groups of male rats (325-375 grams) were tested in the moming (08:00-10:00) 
or in the aftemoon (15:00-17:00) over 5 consecutive days. 

2. Testing for process long-term stability. Adult groups (4 rats/group) of male 
rats (325-375 grams) were tested for flinching behavior over 8 weeks, with separate 
groups being examined at 3-4 day intervals during this period. Adult groups (4 
rats/group) of male rats (325-375 grams) were tested for flinching behavior at 3 groups 
per week for 4 weeks to evaluate process control and capability ming statistical 
process control analysis. All studies were carried out between 08:00 and 10:00 hours. 

3. Testing for body weight effects on flinch behavior. Groups of adult male 
Holtzman Sprague Dawley rats weighing eitherlOO-125 grams (small), 300-350 grams 
(medivim) or 400-450 grams (large) were tested for flinching behavior. No test, control 
article or vehicle was given. 

4. Test of intrathecal catheterization effect on flinching behavior. Groups of 
rats were prepared with lumbar intrathecal catheters. After a 5-day recovery period, 
unimplanted animals of comparable weights (325-350 grams) were tested for formalin 
response. 

5. Test for intrathecal and intraperitoneal drugs on flinch behavior. Groups of 
rats received intrathecal or intraperitoneal injections of morphine sulfate or MK80L 
Dose response curves of each agent were generated to relate drug effect to paw flinch 
and for comparison with previously reported data. 

6. Test of correlation between computer and human observer in detecting 
flinch behavior A trained observer counted flinches each minute for 60 minutes at the 
same time as the automated device performed its data logging. A rat was done once in 
the moming for 6 days. 
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Example 1 

Formalin induced behavioral response 

The injection of formalin into the dorsum of one hind paw of 100 male 
Sprague Dawly rats over a period of 10 months showed a reliable biphasic flinching 
5 of the injected paw, with peak flinch rates during Phase I and II being on the order 
of 50 flinches/min and 37 flinches/min. Figure 7 presents the time course of flinch 
responses, the mean cumulative flinches observed by Phase (Phase I; 203 ± 9 
SEM; H: 1058 ±39 SEM; IIA: 792 ± 27 SEM; and IIB: 266 ± 16 SEM); and the 
statistical distribution of flinching in Phase I and H. Analysis of this distribution 
1 0 (Kolmogorov-Smimov) indicates that each Phase was distributed normally 
(p>0.9999). 

System stability and reliability 

Separate groups of 4 male Holtzman Sprague Dawley rats were run in the 
moming (AM) and aftemoon (PM) on 5 consecutive days. The following 
1 5 summarizes comparisons of the data. 

Comparison between test chambers 

There was no difference in the scores for Phase 1, 11, IIA or DB generated in 
each of the 4 chambers when examined over the 5 days of AM and PM testing 
(Table 2). There were no differences across test chamber for any phase (1-Way 
20 ANOVA,p>0.1 toO.97). 

Comparison of flinching behaviors measured over 5 day interval 

Examination of flinching behavior over 5 consecutive days by phase 
revealed no systematic differences for Phase I, II, HA or IIB (Table 3). There were 
no differences across days for any phase (1-way ANOVA, p>0.07 to 0,32). 

25 
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Table 2. Cumulative flinehing behavior (Mean ±SD) by test chamber, by phase. 















Test Chamber: 












1 


10 


93 ±47 


794 ± 332 


584±215 


210 ±156 


2 


10 


140 ± 35 


740 ±192 


548 ±134 


192 ± 91 


3 


10 


151 ±64 


692 ±331 


508 ± 208 


184 ±149 


4 


10 


131 ±65 


737 ± 162 


542 ±107 


195 ± 79 



Paw flinch system stability measurements 

Table 4 shows the results of the initial data analysis taken on selected days using 4 
animals per day over a period of 7 weeks. For botii Phase I and Phase n the following 
values were calculated: a) the Xbar chart Center Line value, an estimate of the process 
average; b) the Xbar chart Control Limits, which identify the distribution range within 
which the process; c) the S chart Center Line Value, an estimate of the process variation; 
d) the S chart Control Limits, which identify the distribution range within which process 
variation can be considered in "control"; e) the % Upper and Lower Specification Limits 
(%U/L SL) indicating the nmnber of times the measured values exceeded either upper or 
lower control limits and displayed as a % of the total number of tests run; and f) the 

Table 3. Cumulative flinching (Mean ±SD) observed over 5 consecutive days, 
by phase. 















Day 1 


8 


133 ±56 


660 ±175 


522 ±106 


138 ±106 


Day 2 


8 


103 ±59 


620 ±187 


449 ± 93 


171±111 


Day 3 


8 


150 ±60 


889 ±401 


658 ±234 


231 ± 169 


Day 4 


8 


149 ±61 


860 ±175 


609 ±139 


251 ± 73 


Day 5 


8 


109 ±42 


676 ±208 


491 ± 171 


186±110 
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Capability Index (Cp), an indication of whether the formalin testing process was in 
"control" and capable" of staying within the limits specified. 



Table 4. Paw flinch system stability summary. 




Phase I 



257 



±96 



59 



±74 



456 



50 



0.0 



0.0 



1.064 



Phase n 



1298 



±447 



274 



±348 



2186 



353 



0.0 



0.0 



1.026 
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Example 2 

Characterization of the variables influencing the formalin response 

Morning (AM) versus afternoon (PM) test periods. 

Examination of flinching behavior in the morning (08:00-10:00) and in the 
afternoon (15:00-17:00) over 5 days of testing showed no difference between AM 
and PM for all four phases (Table 5). There were no differences across time of day 
for any phase (1-way ANOVA, p>0. 10 to 0.99). 



Table 5. Cumulative flinching behavior (Mean ± SD) as a function of testing 
time, by phase. 
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Effects of body weight on flinching behavior . 

Three groups of rats differing in body weight (125-175 grams, 300-350 
grams and 400-450 grams) were examined for flinching behavior. No difference in 
flinch response was observed in Phase I, n and IIA (1-way ANOVA, p>0.05). A 
weight-related difference was seen in the flinching behavior for Phase IIB (1-way 
ANOVA, p<0.05) (Table 6). 



Table 6. Cumulative flinching behavior (Mean ± SD) as a function of body weight, by 
ihase. 




125-175 grams 



193 ±45 



768 ± 208 



686 ±174 



81 ±69 



300-350 grams 



133 ±56 



660 ±175 



522 ±106 



138 ±106 



400-450 grams 



8 



232 ± 60 



905 ±196 



689 ± 162 



216 ±96 



p >0.05 and * p<0.05, 1-way ANOVA 



10 Effects of gender on flinching behavior. 

Gender differences in flinching behavior was tested using 125-150 gram 
Holtzman Sprague Dawley rats. No difference in flinch response was observed in 
Phase I, n and IIA (1-way ANOVA, p>0.05). A gender related difference was 
seen in the flinching behavior for Phase IIB (1-way ANOVA, p<0.05). (Table 6.) 

15 Table 7. Cumulative flinching behavior (Mean ± SD) as a function of gender, 
by phase. 















Female 


7 


161 ±83 


865 ± 262 


675 ±168 


191±111 


Male 


8 


193 ±45 


768 ±208 


686 ±174 


81 ±69 



p >0.05 and * p<0.05, 1-way ANOVA 
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Table 8. Cumulative flinching behavior (Mean ± SD) as a function of injection 
dose, by phase. 













No Injection 


8 


92 ±56 


146± 102 


92 ±48 


54 ±79 


Saline 


12 


109 ±75 


302± 175 


201 ± 123 


102 ±73 


Formalin, 0.5% 


12 


127 ±66 


400± 185 


301 ± 169 


100 ±79 


Formalin, 1.0% 


12 


171 ±55 


61 3± 265 


495 ±211 


118 ±100 


Formalin, 2.5% 


12 


214 ±75 


898±266 


662 ±210 


236 ±148 


Formalin, 5.0% 


12 


230 ± 84 


81 1± 194 


659 ± 131 


152 ±88 



Effect of formalin concentrations. 
5 The injection of saline into the paw resulted in a modest incidence of 

flinching behavior while formalin concentrations from 0,5%/50^L to 5.0%/50^L 
resulted in a concentration-related increase in flinching behavior measured in all 
phases and as compared to saline (Table 8) (1-way ANOVA, p<0.002). Post hoc 
analysis between groups disclosed phase and concentration differences in flinch 
1 0 response (Table 9). 

Effects of intrathecal catheterization on flinching behavior 

No statistically significant differences were observed in flinch response 
between animal groups with and without intrathecal catheters over all 4 study 
phases (1-way ANOVA, p>0.2). (See Table 10) 

15 Comparison between computer and human observer 

On different days 6 male Holtzman Sprague Dawley rats (300-350 grams) were 
injected with 50 \iL of 5% formalin and an experienced human observer counted 
flinching behavior concurrent with signal acquisition and flinch detection by 
computer. Phase I counts were observed to be lower for the human observer but 
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Table 9. Post hoc group comparisons of formalia injection concentrations, by 
ihase. 




No Injection vs. Saline 



No Injection vs. Formalin, 0.5% 



No Injection vs. Formalin, 1.0% 



No Injection vs. Formalin, 2.5% 



No Injection vs. Formalin, 5.0% 



Saline vs. Formalin, 0.5% 



Saline vs. Formalin, 1.0% 



Saline vs. Formalin, 2.5% 



Saline vs. Formalin, 5.0% 



Formalin, 0.5% vs. Formalin, 1.0% 



Formalin, 0.5% vs. Formalin, 2.5% 



Formalin, 0.5% vs. Formalin, 5.0% 



Formalin, 1.0% vs. Formalin, 2.5% 



Formalin, 1 ,0% vs. Formalin, 5.0% 



Formalin, 2.5% vs. Formalin,5.0% 



S = Significant difference between groups, p<0.05, Tukey/Kramer post hoc test. 



Table 10, Cumulative flinching behavior (Mean ± SD) as a function of 
catheter, by phase. 




total counts remained close between the two methods (Figure 9). A scattergram 

2 

plot of simultaneously acquired data revealed good correlation (r =0.94). 
Example 3 

Effects of systemic and spinal drugs on formalin evoked flinching behavior 

5 Morphine . 

Intrathecal and systemic injection of morphine {\i opioid agonist) resulted in 
a potent dose dependent reduction of all phases of the flinching response 
(FigurelO). The ED50 values were calculated using least squares linear regression 
analysis of the dose effect (% Maximum Possible Effect, %MPE) curves (Table 
10 11). 

MK80L Intrathecal and systemic injection of MK801 (a noncompetitive NMDA 
receptor antagonist) resulted in a potent dose dependent reduction in Phase 11 but 
not Phase I of the flinching response (Figurell). The ED50 values were calculated 
using a least squares linear regression analysis of the dose effect (% Maximum 
1 5 Possible Effect, %MPE) curves (Table 11). 



Table 11. ED50 values for the effects of systemic morphine, intrathecal 
morphine and MK801 on flinching behavior (Mean and 95% confidence 
intervals)'*' 









IT Morphine (^g) 


1.6 (0.6-4.7) 


1.5(0.6-3.6) 


IP Morphine (mg) 


4.7 (7888-.003) 


1.8(1.1-2.9) 


IT MK801 (^ig) 


9.5(87-1.0) 


7.3(43-1.2) 



* Data for IP MK801 were not calculated from dose effect curves due to 
20 lack of monotonic sequence. 
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Power analysis 

The difference in flinch counts necessary to show statistical significance 
was determined using the following assumptions: a) 2 groups being compared with 
8 animals per group; b) a two-tailed analysis; c) a 0.05 level of rejection; and d) a 
reasonably high probability (70% to 90%) of detecting true group mean 
differences. Determinations were based on mean and standard deviation averages 
were derived from 5 different 8-animal formalin control studies (Table 12). 



Example 4 

Injection of formalin into the paw leads to a biphasic flinching behavior with 
10 the magnitude of the behavior positively co varied with the concentration of 
formalin. 

Table 12. Test Group mean values (X) with the power* (70%-90% 
probability) of showing statistical significance from Control means, by phase, 
and expressed as a percent of the Control mean and as a range of applicable 
1 5 values. N=8, p=0.05, two-tailed. 




Phase I 



56 



42%, 
75>X>183 



47%, 
69>X>189 



52%, 
62>X>196 



Phase n 



741 



229 



30%, 
521>X>961 



33%, 
496>X>986 



37%, 
466>X>1016 



Phase HA 



546 



149 



26%, 
403>X>689 



29%, 
387>X>705 



33%, 
367>X>725 



Phase DB 



195 



114 



56%, 
85>X>305 



63%, 
73>X>317 



70%, 
58>X>332 



Power is 1- probability of committing a Type II error. 
'Minimum detectable difference analysis from Biostatistical Analysis. Zar, 1984. 
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Formalin-evoked flinching: Supraspinally organized complex behavior 

From the perspective of a model of nociceptive transmission, an important 
question is whether the formalin-evoked flinching behavior reflects an endpoint that 
is mediated through the exaggeration of spinal traffic to supraspinal centers, or does 
5 the biphasic (Phases I and 11) flinching behavior represent an exaggerated spinal 
reflex? Several observations are relevant. First, injection of formalin induces a 
variety of complex, unconditioned behaviors, such as licking and guarding of the 
injected paw, which are reflective of a higher order motor organization. Second, 
spinal transections markedly reduce Phase I flexion-extension of the formalin 

10 injected hind paw and Phase n is abolished (42). Third, in addition to the paw 
withdrawal, formalin injection induces biphasic autonomic (cardiovascular) and 
supraspinally mediated hormonal responses (29, 38, 51). Fourth, assessment of the 
firing pattems of dorsal horn wide dynamic range neurons, many of which are 
projection neurons, has shown that comparable formalin injections results in a 

1 5 similar biphasic activation pattem ( 1 5, 33) and an increase in the expression of 
cFOS in dorsal horn neurons (1). Jointly, these findings support the argument that 
Phase n behaviors reflect augmented responses to spinofligal traffic and not simply 
an augmented reflex. 

Role of small afferents in formalin evoked flinching. 

20 The flinch response is evoked and maintained by persistent small afferent 

input. Two lines of evidence support this hypothesized mechanism: i) treatment 
with the C-fiber neurotoxin capsaicin reduces the response to an irritant injected 
into the paw, suggesting a role for small afferents (17, 29, 46); and ii) local 
anesthetic blockade of the afferent input during Phase II reduces dorsal hom neuron 

25 activity (15) and halts flinching and grooming (13, 39) during Phase I and Phase 11. 
In other words, all flinching behavior (Phase I and Phase 11) requires ongoing 
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afferent traffic. 

Mechanisms underlying the biphasic components of the flinching behavior 

Measurement of the firing pattern of small sensory afferents evoked by 
formalin injected into their receptive fields reveals an acute burst of activity that 
5 persists for several minutes. This initial discharge is followed by a persistent low 
level of afferent activity in slow (small diameter) and fast (large diameter) sensory 
afferents (32), Temporally, the initial intense flinching (Phase I; 0-10 min) 
correlates with the initial afferent barrage, while the prominent second phase of 
flinching (Phase 11: 20-60 min) corresponds with the interval when there is a 
10 relatively modest, but nonzero level, of afferent input. An important question is, 

given that there is reduced afferent traffic during Phase 11 (32), what is the origin of 
the prominent flinching that is observed during Phase II? We believe that the 
system is unexpectedly complex, involving several potential mechanisms, all of 
which may contribute to this second phase. 

15 Peripheral components . 

As noted above, based on local anesthetic blockade, ongoing afferent traffic 
is essential for dorsal hom neuron activity and the delayed onset of behavior during 
the second phase afl;er formalin. Peripheral inflammation generated by formalin can 
initiate activity, perhaps in populations of small cutaneous afferents that are not 
20 normally active and have the ability to strongly drive dorsal hom neurons ("silent 
nociceptors") (19). This provides a peripheral mechanism that would initiate 
enhanced neuronal activity and flinching behavior in the face of an apparent 
reduction in overall afferent traffic during the second phase (13). 

Central Components. 

25 Persistent small afferent input evokes a facilitation of spinal nociceptive 
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processing (14, 44). Accordingly, it is reasonable to conclude that intradermal 
formalin may initiate such a cascade, leading to a state of central facilitation that is 
maintained during the low level of afferent traffic in Phase IL Such a facilitated 
state, combined with a low level of ongoing afferent traffic, would provide a 
5 mechanism for the observed high level of flinching during Phase 11 (14). Two 
observations support this role of a central sensitization: i) delivery of agents 
believed to suppress small afferent input (e.g., opiates) reduces the Phase II 
response, even when their action is limited to the interval of Phase I (3, 9), but see 
studies with systemic opiates (40); and ii) delivery of classes of agents that do not 
1 0 block acute excitation or acute pain behavior, but block afferent-evoked spinal 
facilitation such as NMDA antagonists, COX inhibitor and NOS inhibitors, 
diminish the second phase of the formalin-evoked behavior (see Yaksh, 1997 (45) 
for references). 

Interactive contributions. 

15 It is believed titiat it is likely that changes in both afferent traffic (peripheral) 

and spinal processing (central sensitization) components contribute to the observed 
formalin-induced behavioral states. As noted, sensory afferent recordings have 
emphasized that formalin injection leads to ongoing activity in large and small 
afferents. Large afferent activity (e.g. light touch) is not typically associated with 

20 the initiation of a pain state. However, after local injury or small afferent 

activation, an exaggerated response of dorsal horn neurons to high intensity stimuli 
applied to the site and to low intensity tactile input applied adjacent to the injury 
site is detected (see Dickenson et al. (14)). These phenomena have the behavioral 
parallel of 1 ° hyperalgesia and 2° tactile allodynia (8, 28). Pharmacological studies 

25 have shown that the 1 hyperalgesia has a well defined peripheral component, while 
the 2' tactile allodynia is initiated but not sustained by the primary injury input. We 
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thus hypothesize that the 2nd phase pain behavior may arise from: i) an ongoing 
central sensitization initiated and maintained by small afferent input; and ii) the 
afferent traffic coming in from small afferent (normally high threshold) nociceptors) 
and perhaps of equal importance from large afferents (low threshold, 
5 mechanoreceptors) that can induce pain behavior when there is a central 
sensitization. 

Manual Formalin Testing 

The evident utility of the formalin response suggests that it serves as a robust 
model for screening anti-hyperalgesic agents. A limitation is the time and training 

1 0 required to perform the test. In previously published work, the primary behavioral 
index has been the counting of the number of flinches in which 1 minute flinching 
epochs are counted every 5 or 10 minutes over a 60 minute observation interval. 
This spacing of observations permits a single technician to follow several rats 
concurrently. Other methods, such as the counting of time spent in multiple 

1 5 behaviors ( 1 2, 42), might provide additional sensitivity, however, the complexity of 
such an analysis places limitations on the implementation of the model: i) It reduces 
in the number of animals that can be followed concurrently; ii) the need for 
observer vigilance means that with repeated testing during the workday, the drift in 
operator reliability can surely be anticipated. It is important to note that in work 

20 published on the method, such practical issues as reliability between observers over 
time are not discussed; and iii) the use of indices that involve judgments on weight 
bearing or appropriate licking behavior requires that the observer undergo extensive 
training and some form of technique validation to prove that needed skills have 
been acquired before considering that the data produced is useable. Again, most 

25 investigators do not generally discuss these issues. Wheeler-Aceto and Cowan (42) 
specifically cited the simple counting of formalin flinches as reducing the degree of 
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inter-observer variation. Coderre et al. (12) emphasized the degree of inter-observer 
reliability but, again, this is Ukely a best case analysis unencumbered by drug or 
behavior related perturbation. 

Automated systems 

5 Given the limitations associated with manual assessment of formalin-evoked 

flinching, efforts have been made to develop automated assessment systems. The 
systems described include one by Jett and Michaels that involves measurement of 
mass shifts secondary to the movement of a rat in a confined cylindrical cage (21). 
The use of a complex filtering algorithm is said to reduce the contribution of v^hole 

10 body movement from those arising from the movement of the much smaller mass of 
the hind paw. This model has been employed in several reports (7, 29). Another 
model utilizes a video camera system and relies on pattern recognition. The 
algorithm employed is not described, but it appears to depend upon the symmetry of 
the animal (22). In neither case can the movement of the injected paw be defined 

1 5 with certainty. In contrast, the model presented here addresses this issue directly. 

Modeling of the detector system 

At the outset, the aim was to minimize radio frequency interference and the 
interaction of the coil fields with the surrounding environment. Collar size was 
configured to be unencumbering to the rat and constructed of a metal with 

20 sufficient permitivity to generate a useful signal. Modeling of the electromagnetic 
field and collar interaction was undertaken to determine a workable system 
configuration and a configuration was defined that uses an open collar to produce a 
usefiil signal within a relatively localized field at the frequency range of 6-8 KHz. 
Computational analysis indicates that increased field strengths can be achieved by 

25 using higher coil voltage, or by using a fixed signal strength, and: i) metals of 
higher permitivity; ii) a higher sinusoidal frequency; and/or iii) increasing field 
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density by decreasing the field coil size. 



Validation of the automated formalin system 

The two methods used here to determine system relevance were: obtaining 
human and system flinch counts concurrently, and comparing flinch count data 
5 obtained from previous drug studies with similar new studies using the automated 
system. 

As indicated, concurrent paw flinch assessment by a trained human observer 
and by the automated system demonstrated good correlation. There was a tendency 
for the automatic flinch counts in Phase I to be higher than those obtained by the 
1 0 observer, but this is believed to reflect the inability of the human observer to follow 
the rapid flinching which is evident during that phase. 

With regard to the activity of drugs known to alter flinching behavior, 
morphine given intrathecally and systemically produced a dose-dependent reduction 
in all phases of the response to formalin. These results are consistent with the 

1 5 mechanism of action of morphine in which at the spinal level, \i opiate receptors are 
believed to diminish the excitability of input associated with the activation of small 
primary afferents (45). The ED50s obtained from the automated assessment 
correspond closely with those ED50s reported using human scoring. Intrathecal 
delivery of an NMDA receptor blocker, MK801, has been shown to block Phase 11 

20 of the response to formalin and, to a lesser degree. Phase I in a dose- dependent 
manner. Interestingly, while the intrathecal response curve was clearly monotonic, 
there was an apparent increase in flinching behavior at the lowest doses after 
systemic delivery of the drug. Examination of the animals during the study indeed 
confirmed that there was enhanced activity leading to increases in flinch count. 

25 The reason for this biphasic component is not known, but it is appreciated that the 

non-competitive NMDA antagonists may have stimulatory properties that might 
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reflect a PCP-like action. Such effects may be accounted for if the systemically 
delivered drug is getting into the brain (as compared to a spinal action). 



Assessment of system reliability 

In any assay system, base line stability over time is a primary concern. 
5 Baseline flinching behavior was assessed in separate studies over a 5 day, 7 week, 
and 1 1 month time increments. Over these intervals there were no statistically 
significant shifts in baseline activity in standard untreated animals. 

A more definite approach is to undertake a statistical process analysis. 
Using a statistical process control analysis, we specijBcally focused on the 7 week 

0 period. We sought to determine if over the 7 week time increment the process, 
measurement of paw flinches, was in control (i.e., paw flinch measurement varied 
only within the limits of a selected statistical distribution) and capable (i.e., a high 
percentage of the measurements taken fell within previously determined 
specification limits). A process that is "in control", or stable, can be defined by a 

5 normal distribution of the observations obtained over time. One measure of 

process control is the capability index (Cp). This measure defines how satisfactory 
a process is at meeting the requirements placed on it. A Cp value of greater than 
1 .00 means that only a small number of the formalin tests would not fall within 
the limits considered relevant to the normal testing process. Conversely, a Cp<LO 

0 indicates that a significant number of observations fell outside of the specified 

distribution. In the present studies, Cp values for both Phase I and II over a 7 week 
interval were > 1 .000 and thus "in control" and "capable". 

Though seldom used in behavioral system studies, these types of analyses 
can: i) bolster confidence in a testing process, ii) help to define its sensitivity 
'5 limits, and iii) aid in the detection of tr^sient or systematic changes before they 
affect study results. 
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Example 5 

Variables influence the formalin response 

Several variables may impact upon the magnitude of the flinching behavior. 
Using the automated assessment systems^ selected variables were examined, 

5 Behavioral indices . 

Injection of formalin into the paw results in a variety of spontaneous 
behaviors including flinching of the paw, licking of the paw, lifting the paw from 
the surface of the chamber, changes in weight bearing, and vocalization. 
Assessment of the intensity of a behavior induced by formalin injection has 

10 typically involved one or a combination of measures including the measurement of 
duration and/or frequency for one or more of the known behaviors (2, 12, 1 8, 42). 
In most cases, as formalin concentration increases or as the dose of an analgesic 
agent (e.g., morphine) decreases, there is a corresponding increase in the assessed 
index (see references in Table 13), Several authors have demonstrated that 

15 weighted measures (e.g., time of licking and flicking) may provide a more robust 
index with a greater dynamic response at lower formalin concentrations than other 
simple indices such as flinching. Choice of indices may be influenced by the time 
required to complete a test and the level of training required to achieve reliability. 

Concentration of formalin. 

20 The current literature indicates that increasing formalin concentration, 

typically over a range of 0.5 to 5%, leads to a more intense and progressively longer 
lasting flinch response. Assessment of the intensity of several pain behaviors have 
shown that as formalin concentration rises, the magnitude or incidence of the 
measured behaviors, blood pressure or flinching increases (2, 5, 1 1, 12, 39, 41). In 

25 the present studies, we observed a plateau effect such that the maximum was 
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observed at 2.5%. These results, obtained with automated flinch counting, are in 
accord with previous work using weighted behavioral measures and blood pressure 
(11.12.39). 

Rat strain, 

5 The present studies were done with Sprague Dawly, Holtzman strains. There 

are few systematic studies describing strain differences with regard to formalin 
response. A variety of rat strains have been used including Sprague Dawley, Harlan 
strain (25) and Sprague Dawley, Holtzman strain (see Dirig and Yaksh (16); Wistar 
(2) and Long Evans (12)). Taylor and colleagues noted that the magnitude of the 
1 0 blood pressure response observed during Phase I was similar in Sprague Dawly rats 
obtained from two suppliers (Charles Rivers and Bantin Kingman), but that the 
latter displayed a significantly reduced Phase II response. 

Bodv weight / Age/Gender 

Body weight varies directly with age. In the present work, young (1 10-125 
1 5 grams) rats flinched less tiian larger, older, rats (400-450 grams) of the same strain. 
There appear to be little differences in flinching behavior as a fimction of gender. 
In our work, when matched for size, females and males of a single strain showed 
comparable flinching behaviors, with the exception of the late component (Phase 
IIB) in which females showed a statistically greater response than males. The 
20 significance of this modest difference is not known. 

Catheterization. The surgical placement of a spinal catheter had no 
significant effect upon flinching behavior as compared to the unimplanted rat. 
been examined for their effects upon Phase 11 of the formalin test after systemic and 
intrathecal delivery. Agents may be functionally considered in two classes, those 

25 
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Example 6 

Pharmacology of the Formalin Test 

The pharmacology of the formalin test has been a subject of considerable 
investigation. Table 13 summarizes a number of families of agents, which have 



Table 13. Effects of intrathecal and systemically delivered pharmacological 
of a|!ent^ on '''^^I^^^ ^^^^ El^^^^^. 



Agonists 






\i opiate 


(25) 


(42) 


d opiate 


(25) 




K opiate 


(25) 


(42) 


Alpha 2 adrenergic 


(25) 




Adenosine A- 1 


(25) 




Nicotinic 




(6) 


GABAA 


(16) 




GABA-B 


(16) 




Benzodiazepine 


(16) 




Gabapentin 


(34) 


(36) 


Antagonists 






NMDA antagonist 


(10) 


(20) 


Glycine site antag 




(20) 


Kainate antagonist 




(35) 


AMPA antagonist 


(10) 


(20) 


NK-1 antagonist 


(50) 


(21) 


COX inhib 


(24) 


(24) 


EP antagonists 


(25) 




NOS inhibitor 


(26) 




Sodium ch blocker 




(4) 


N-Ca ch blocker 


(27) 




Cholinesterase 
inihib 


(30) 





37 



that will completely reduce Phase 11 (the facilitated component) and those that 
appear to exert a significant, but limited, plateau effect. To date, all agents that are 
effective for acute pain processing fall into the first category and include opiate and 
5 alpha 2 agents. Agents which fall in the second category include, for example, 
NMDA and NKl receptor antagonists, Adenosine Al agonists, GABA A/B 
agonists, and cyclooxygenase inhibitors (47, 48). hi this invention, we have 
demonstrated that the automatic assessment of flinching behavior provides data that 
are both qualitatively and quantitatively similar to data that has been previously 
1 0 reported by human observers. 

CONCLUSIONS 

The current literature suggests that the formalin test reflect a neural substrate 
that involves the generation and support of a facilitated state of processing. This 
behavior observed during Phase n appears to arise out of the initial barrage of 
15 afferent traffic and continuous low-level input found in a formalin-induced injury. 

The response of this system is such that it appears to reveal a pharmacology 
that is considered to reflect facilitated states of processing. 

An automated sensing system has been developed that counts the spatial 
displacement of the injected hind paw. This system has been shown to describe a 
20 response count and distribution that resembles that obtained with manual counting 
systems. 

The limitations of this approach are that it only counts paw flinching and 
does not assess time spent in other behaviors deemed to be representative of the 
animals pain behavior. Nevertheless, based on sensitivity to formalin 
25 concentrations, duration of behavioral responses and response to drugs, there 
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appear to be no distinguishable differences between simple flinch analysis and a 
weighted behavioral score. 

The strength of the system lies in the ability to screen large numbers of 
animals on a daily basis, absence of operator fatigue, a likely increase in reliability 
5 over repeated tests and reduced training time. 

In accordance with the preceding explanation, variations and adaptations of 
the dipolar device for presensing impending cardiac disturbances of the present 
invention will suggest themselves to a practitioner of the medical instrumentation 
arts. In the spirit of this invention, these and other possible variations and 
10 adaptations of the present invention, and the scope of the invention, should be 
determined in accordance with the following claims, only, and not solely in 
accordance with that embodiment within which the invention has been taught. 
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